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ABSTRACT 

By means of self-consistent three-dimensional (3D) magnetohydrodynamics (MHD) numerical simu- 
lations, we analyze magnetized solar-like stellar winds and their dependence on the plasma-/? parame- 
ter (the ratio between thermal and magnetic energy densities). This is the first study to perform such 
analysis solving the fully ideal 3D MHD equations. We adopt in our simulations a heating parameter 
described by 7, which is responsible for the thermal acceleration of the wind. We analyze winds with 
polar magnetic field intensities ranging from 1 to 20 G. We show that the wind structure presents 
characteristics that are similar to the solar coronal wind. The steady-state magnetic field topology for 
all cases is similar, presenting a configuration of helmet streamer-type^ with zones of closed field lines 
and open field lines coexisting. Higher magnetic field intensities lead to faster and hotter winds. For 
the maximum magnetic intensity simulated of 20 G and solar coronal base density, the wind velocity 
reaches values of ~ 1000 km s~^ at r ~ 20 ro and a maximum temperature of ~ 6 x 10^ K at r ~ 6 tq. 
The increase of the field intensity generates a larger "dead zone" in the wind, i. e., the closed loops 
that inhibit matter to escape from latitudes lower than ~ 45^ extend farther away from the star. The 
Lorentz force leads naturally to a latitude-dependent wind. We show that by increasing the density 
and maintaining Bq = 20 G, the system recover back to slower and cooler winds. For a fixed 7, 
we show that the key parameter in determining the wind velocity profile is the /^-parameter at the 
coronal base. Therefore, there is a group of magnetized fiows that would present the same terminal 
velocity despite of its thermal and magnetic energy densities, as long as the plasma-/? parameter is 
the same. This degeneracy, however, can be removed if we compare other physical parameters of the 
wind, such as the mass-loss rate. We analyze the infiuence of 7 in our results and we show that it is 
also important in determining the wind structure. 

Subject headings: stars: winds, outfiows - MHD - methods: numerical - stars: late-type 



1. INTRODUCTION 

Studies of the solar corona (SC) have played a crucial 
role in understanding stellar winds in general. Due to 
our privileged position immersed in the solar wind, re- 
searchers have had access to a great quantity of data that 
allow a detailed understanding of the physics that is op- 
erating in the Sun. Recent sophisticated observations, e. 
g., via SOHO and Ulysses (McComas et al. 1995; Suess 
& Smith 1996; Jones et al. 1998; Wilhelm 2006, among 
others), have shown that the SC is a highly complex sys- 
tem. It consists of long-lived features, like the fast and 
slow solar wind, streamers and coronal holes, and also 
of short-lived features, like the coronal mass ejections, 
solar fiares and sun-spots. Hence, if one could separate 
the long-lived structures from the short-lived ones, one 
could, in principle, come to a better understanding of 

Electronic address: ahne@astro.iag.usp.br 



the SC itself and other similar outfiows. In the present 
paper, we aim to investigate in detail the behavior of the 
long-lived coronal features in solar-like stellar winds. 

Direct measurements of tenuous coronal winds for 
other stars rather than the Sun have proved to be very 
difficult to do, although indirect detections of stellar 
coronal winds have been performed (Wood et al. 2005). 
As it occurs in the SC, it is very probable that the mag- 
netic field is playing an important role in coronal winds 
of solar-like stars. Magnetic activity has not only been 
detected in other stars (Bagnulo et al. 2002; Reiners & 
Basri 2006; Phan-Bao et al. 2006; Wade et al. 2007), but 
also in stars similar to the Sun (Robinson et al. 1980; 
Petit et al. 2005). 

In the absence of a magnetic field, a non-rotating 
corona expands spherically (Parker 1958; Velli 1994). A 
different picture, however, is expected if the magnetic 
energy density at the base of the corona is at least of the 
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same amount of the thermal energy density. 

A particular case of equal magnetic and thermal en- 
ergy densities was explored by Pneuman & Kopp (1971, 
PK71 from now on) for the case of the magnetized so- 
lar wind. By means of a numerical iterative method, 
they concluded that the SC is composed of two differ- 
ent magnetic structures in large scale: a region of closed 
magnetic field loops near the star at low latitudes and at 
high latitudes, open field lines that cannot restrain the 
expanding gas. 

Astrophysical outflows have long been studied (Parker 
1958; Weber & Davis 1967; Mestel 1968; Nerney & Suess 
1975). Several analytical studies were made toward the 
understanding of an expanding magnetized corona. For 
instance. Low & Tsinganos (1986) found a class of radial 
analytic solutions for magnetized steady-state winds by 
removing the latitude dependence of the problem. By 
varying the global magnetic field geometry and other 
parameters of their model, such as the velocity at the 
base of the wind, they analyzed the behavior of flows 
along open magnetic field lines. Tsinganos & Trussoni 
(1991) relaxed the assumption of a poly tropic equation 
of state, but neglected the meridional component of the 
flow. Sauty & Tsinganos (1994) found a new class of 
analytical solutions calculating the exact shape of a field 
line along the flow rather than assuming a global mag- 
netic field geometry. They showed that for a rotating 
young stellar object the shape of the field lines is an 
important parameter for the formation of a coUimated 
jet/non-collimated wind. In a more recent work, Lima et 
al. (2001) constructed exact solutions for a rotating, mag- 
netized wind. They verified that the wind is highly non- 
spherically symmetric. However, as they pointed out, 
their model has limitations, as, for instance, meridional 
flows and the meridional component of the magnetic field 
were neglected. 

Meanwhile, numerical studies were carried out with 
increasing level of sophistication. Washimi & Shibata 
(1993) studied the influence of the stellar rotation on the 
wind structure and acceleration performing 2D simula- 
tions. Solar wind parameters were used, except for the 
uniform angular velocity at the surface of the star that 
was varied as to analyze the effects of the centrifugal 
force. They showed that if the angular velocity is in- 
creased more than ten times the solar value, the centrifu- 
gal force becomes comparable to the thermal force, in- 
fluencing the meridional structure of the wind. Keppens 
& Goedbloed (1999) presented detailed one-dimensional 
(ID) and two-dimensional (2D) ideal MHD numerical 
simulations of a polytropic, axisymmetric wind. In a pos- 
terior work, Keppens & Goedbloed (2000) modeled stel- 
lar axisymmetric rotating outflows by solving the ideal 
MHD equations, investigating the effects of open and 
closed magnetic field lines in the wind. By varying the 
extension of the closed field line region, as well as the in- 
tensity of the magnetic field, they showed that the global 
wind structure is modified. Ud-Doula (2003), working 
towards the modeling of hot-star winds, presented a sim- 
ulation of the solar wind using time- dependent, axisym- 
metric, MHD simulation. He showed that the MHD mod- 
eling was consistent with the work of PK71. Modeling 
the propagation of a coronal mass ejection from the in- 
ner solar corona to 1 AU, Groth et al. (2000) numerically 
reproduced the steady-state bi-modal nature of the so- 



lar wind with a prescribed ad-hoc heating mechanism. 
Working on the modeling of the solar wind, Roussev et 
al. (2003) simulated the 3D structure of the solar wind 
under steady-state conditions, using solar magnetogram 
data as input parameters for the initiation of the wind, 
and considering a variable heating mechanism. Cohen et 
al. (2007) extended Roussev et al. (2003) 's work consid- 
ering, as a heating mechanism, a radial dependence of 
the ratio of specific heats, 7 = 7(r), as to reproduce the 
observed bi-modality of the velocity of the solar wind. 

Despite all the notable evolution of both analytical 
and numerical studies performed in the last decades, 
we are far from a satisfactory 3D MHD description of 
a magnetized wind. Several approximations were made 
in order to make the system analytically and numerically 
tractable (e. g., neglecting meridional flows, assuming a 
polytropic equation of state, assuming a magnetic fleld 
topology). 

In the present study, we investigate the influence of 
the magnetic field in solar-like stellar winds with differ- 
ent plasma- /5 (the ratio between thermal and magnetic 
energy densities). We solve the fully 3D MHD equa- 
tions with the temporal evolution of the energy equation. 
Therefore, the topology of the field is not restricted and 
the steady-state arises from the dynamical interplay of 
the outfiow and the field. Also, meridional fiows arise 
naturally in the system. We neglect the stellar rotation. 
The results presented here are thus valid for non-rotators 
or slow rotators in the region where the toroidal compo- 
nent of the field is still much smaller than the poloidal 
component. 

The paper is organized as follows. In §2, we present the 
numerical scheme used and in §3, the results obtained. 
§4 is dedicated to conclusions and discussion. 

2. THE NUMERICAL MODEL 

To perform the simulations, we make use of the Block 
Adaptive Tree Solar-wind Roe Upwind Scheme (BATS- 
R-US), a three-dimensional MHD numerical code devel- 
oped at the Center for Space Environment Modeling at 
University of Michigan (Powell et al. 1999). 

BATS-R-US uses a computational domain that is 
block-based, consisting of Cartesian blocks of cells that 
can be adaptively refined. It has been widely used to sim- 
ulate the Earth's magnetosphere (Ridley et al. 2006), the 
heliosphere (Roussev et al. 2003; Cohen et al. 2007), the 
outer-heliosphere (Linde et al. 1998; Opher et al. 2003, 
2004), coronal mass ejections (Manchester et al. 2004; 
Lugaz et al. 2005), and the magnetosphere of planets 
(Toth et al. 2004; Hansen et al. 2005), among others. In 
this work, we adapted the version for the outer helio- 
sphere (Opher et al. 2003, 2004) to study the problem of 
the wind of a solar-like star. 

The code solves the ideal MHD equations, that in the 
conservative form are given by 



dt 



| + V.(.u) = 
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puu-\- [p 
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(1) 



PS (2) 
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de 

dt 



(u-B)B 



= pg-u, (4) 



where p is the raass density, u the plasraa velocity, B the 
magnetic field, p the gas pressure, g the gravitational 
acceleration due to the central body, and s is the total 
energy density given by 



s = 



pu^ 



7 



p_ 

- 1 Stt 



(5) 



We consider ideal gas, so p = pksT / {prrip)^ where ks is 
the Boltzmann constant, T is the temperature, prrip is 
the mean mass of the particle, and 7 is the ratio of the 
specific heats. 

Due to the lack of knowledge of all the detailed pro- 
cesses that take place in a stellar wind, it is difficult 
to estimate all the mechanisms that modify the heat 
content of the wind (e. g., conduction, radiation, me- 
chanical dissipation of energy that is transferred to the 
plasma). MHD waves and turbulence are known to play 
an important role on solar wind acceleration and heat- 
ing (Matthaeus et al. 1999; Cranmer & van Ballegooi- 
jen 2005; Cranmer et al. 2007). For this reason, they 
are often included in stellar wind models as well (e.g., 
Airapetian et al. 2000; Vidotto & Jatenco-Pereira 2006; 
Suzuki 2007). The solar-stellar wind connection is usu- 
ally made by scaling the observed solar wind characteris- 
tics (e.g., the wave flux and spectral slope, solar magne- 
tograms) to other stars (Jatenco-Pereira & Opher 1989; 
Schrijver 2001). However, it is not clear how these char- 
acteristics scale to other solar-like stars, in this first anal- 
ysis we start with a simpler treatment of the mechanisms 
of acceleration of the wind. We adopt, therefore, an ap- 
proach similar to Roussev et al. (2003) who considered 
that 7 is associated with "turbulent" internal degrees of 
freedom, in a way analogous to the Sun, where a signifi- 
cant amount of energy is stored in the form of waves and 
turbulent fluctuations. 

2.1. The Grid Adopted 

For all the simulations done, we adopted the same grid 
resolution. Initially, the simulation domain is reflned in 
five levels. Other five refinement levels that are body- 
focused and focused on the equatorial plane (current 
sheet region) are applied next. Finally, an additional 
level is applied to the body. There are 9.1 x 10^ cells 
in the domain. The smallest cell size is 0.018 ro, lo- 
cated around the central body. The maximum cell size is 
4.68 ro- The cell size near the current sheet is 0.036 ro- 
The grid is Cartesian {x^y^z} and the center of the star 
is placed at the origin. The axes x, y and z extend from 
—75 ro to 75 tq. The grid can be seen in Fig. 1. 

2.2. The Initial Conditions 

The simulation domain is initialized with a solar-like 
star at its center with 1 Mq and ro = 1 Rq. The grid 
is initialized with a ID hydrodynamical (HD) wind for 
a totally ionized plasma of hydrogen, he initial temper- 
ature profile is set according to the value of 7 adopted. 
This solution is dependent solely on the choice of the base 
temperature of the wind and the only physical possible 
solution is the one that becomes supersonic when passing 
through the critical radius (Parker 1958). Due to conser- 
vation of mass of a steady wind, we obtain the density 



profile from the radial velocity profile Ur{r). Initially, we 
chose a coronal base temperature of Tq = 1.56 x 10^ K 
and a base density of po = 1.544 x 10~^^ g cm~^ (we vary 
po in the second and third sets of simulations presented 
later on). 

The simulations are initialized with a bipolar mag- 
netic field configuration described in spherical coordi- 
nates {r, 0, cp} by 



^ (cos sin 0,0 



(6) 



where Bq is the magnetic field intensity at the poles, r 
is the radial coordinate, is the co-latitude, and (p is 
the azimuthal angle measured in the equatorial plane. 
The system is then evolved in time until steady-state is 
achieved. 

2.3. The Boundary Conditions 

The inner boundary of the system is considered to be 
the base of the wind at r = ro and its conditions are 
dependent on local fiow conditions: plasma can freely 
leave the reservoir (i.e., the base of the coronal wind), 
but no "backflow" is allowed. Fixed boundary conditions 
were adopted at r = ro . The outer boundary has outflow 
conditions, i. e., a zero gradient is set to all the primary 
variables. 

3. RESULTS 

Table 1 presents the parameters used for the simula- 
tions performed. SOO is a purely HD simulation included 
for comparison purposes. The simulations are divided in 
three sets. The first set of simulations, composed by sim- 
ulations SOI to S05, aims to investigate the effect of the 
magnetic energy density on the wind. We increased Bq 
from 1 to 20 G while maintaining the other initial values 
fixed: Po, the ratio between the surface thermal pressure 
Po and the magnetic pressure evaluated at the pole 



Bo' 



(7) 



is decreased from 1 to 1/400. SOI is similar to the case 
studied by PK71 with /^o = 1. The initial HD wind is 
the same for all the simulations of the first set {po = 
Pe = 1.544 X 10-1^ g cm-^ To = 1.56 x 10^ K, as can 
be seen in Table 1). In the second set of simulations, 
from simulations S06 to S09, Bo is fixed at 20 G and 
we vary po from 1 to 400 times the solar density, pc = 
1.544 X 10~^^ g cm~^, thus increasing Po back to 1 in 
S09. In both sets of simulations, we adopted the value 
7 = 1.01. We then consider a third set of simulations, 
where parameters similar to the first and second sets are 
used, but we adopt 7 = 1.1. 

3.1. The First Set of Simulations 
3.1.1. Evolution to Steady- State 

The Lorentz force Fb acting on the ffow is given by 
the difference between the magnetic tension, 7b, and a 
non-isotropic gradient of the magnetic pressure Pb 

Fb = Tb - VxPb , (8) 

where V_l is the component of the gradient perpendicu- 
lar to the magnetic field B. As the magnetic fleld lines 
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Fig. 1. — Three different views of the grid. Inside each region, we indicate the maximum cell size of the region (in ro units). 



TABLE 1 
The complete set of simulations 



Name 


To (MK) 


Po/Pc* Bo (G) 




7 


SOO 


1.56 


1 


oo 


1.01 


First Set of Simulations 


SOI 


1.56 


1 1 


1 


1.01 


S02 


1.56 


1 5 


1/25 


1.01 


SOS 


1.56 


1 10 


1/100 


1.01 


S04 


1.56 


1 15 


1/225 


1.01 


S05 


1.56 


1 20 


1/400 


1.01 


Second Set of Simulations 


S06 


1.56 


1.8 20 


1/225 


1.01 


SOT 


1.56 


4 20 


1/100 


1.01 


SOS 


1.56 


16 20 


1/25 


1.01 


S09 


1.56 


400 20 


1 


1.01 


Third Set of Simulations 


SOlb 


1.56 


1 1 


1 


1.1 


S02b 


1.56 


1 5 


1/25 


1.1 


S08c 


1.56 


16 20 


1/25 


1.1 


S09c 


1.56 


400 20 


1 


1.1 


* 

pc = 


1.544 X 10- 


16 g cni-3 







are more curved at low latitudes, the tension becomes 
more important as one approaches the equatorial plane. 
Hence, opening a bipolar field structure at low latitudes, 
where lines are more curved, is more difficult than open- 
ing a bipolar field topology at high latitudes. As both 7b 
and V_lPb are latitude-dependent, so is Fb- We there- 
fore expect from the interaction of the outflow and the 
magnetic field a latitude-dependent wind. 

Figures 2a and 2b present the temporal evolution of 
the magnetic field lines from iteration N = 1 to N = 
20, 000. This evolution was common to all runs. The plot 
presented is in the meridional plane, for 803. Figure 2a 
shows that at high latitudes, field lines tend to become 
open. Figure 2b shows the initial stretching of the lines 
emerging from low latitudes along the equator. 

As the simulation evolves in time, the magnetic field 
lines reconnect in the equatorial plane. The steady-state 
configuration (Fig. 2c) consists of a formation of a bi- 
component wind. This configuration is composed by a 



zone of open magnetic field lines (I) coexisting with a 
zone of closed loops anchored on the star (II) , also called 
helmet streamer. At the top of the closed loops lies the 
neutral point of cusp- type (III): approaching the neutral 
point from inside the closed loops, the magnetic field goes 
to zero. Beyond the neutral point, along the equatorial 
plane, a current sheet is formed (thick line in Fig. 2c). 
The zone of closed field lines is located at low latitudes 
from ~ —45^ to 45*^ and extends up to the neutral point. 
Beyond the neutral point, the zone of open field lines 
fills all the volume. Figure 3 shows a 3D view of the final 
configuration. 

3.1.2. Steady-State Wind Profile 

Figure 4 shows the total velocity of the wind in the 
meridional plane for the first set of simulations. The 
magnetic field lines are the black lines and the Alfven 
surface is indicated by the white line. It can be seen 
that the increase of Bq leads to faster winds. 

The five panels present similar structures but quanti- 
tatively very different. The wind is not spherical; higher 
velocities are achieved at high latitudes. The higher 
is, the higher is the departure from spherical symmetry 
and the higher is the total velocity of the wind. 

The bi-modality of the wind is due to the nature of 
the magnetic force. A purely HD (non-rotating) wind is 
spherically symmetric but in the MHD case, this symme- 
try is lost because the magnetic force has a meridional 
component. 

The solar wind has an observed bi-modal velocity rang- 
ing from 400 to 800 km s~^ at 1 AU (solar minima). Al- 
though our model does not make use of magnetograms 
and thermal latitudinal heating as more realistic solar 
models such as Cohen et al. (2007) and Liu et al. (2008) 
do, it is interesting to compare SOI to their results. In Liu 
et al. (2008) 's model, at r = 4.5 ro, the highest velocity 
of ^ 440 km s~^ is achieved at the polar axis {0 = 0*^), 
while in the mid-latitudes {0 ~ 45*^), this value is de- 
creased to ~ 230 km s~^. For SOI, the same positions 
leads to 201 km s~^ and 190 km s~^, respectively. In our 
models, we are not treating the presence of waves and 
damping that is responsible for the latitude dependence 
that is observed in the solar wind. However, it is inter- 
esting that, just by the presence of magnetic field with an 
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Fig. 2. — (a) and (b): Evolution of the magnetic field lines for S03. The color of each streamline represents a given iteration, ranging 
from N = 1 {t = 0, blue) to 20, 000 {t = 7h, black). At = 1, the field is a dipole. The steady-state is achieved at = 148, 000 {t = 51h). 
(c): Final steady-state configuration of the magnetic lines: (I) represents the zone of open field lines; (II) the zone of closed field lines; (III) 
the location of the neutral point. The thick line represents the position of the current sheet formed. 



TABLE 2 

Co-latitude dependence of Ur at r = 4.5 ro 

Name nr(6> = 18") nr(6> = 45") Aur^ 

501 200 189 n 

502 243 224 19 

503 329 295 34 

504 456 400 56 

505 610 528 82 




Fig. 3. — 3D view of the steady-state for S03. 



initially spherical wind, we obtain a latitude-dependent 
wind, although not as dramatic as in the Sun. 

Figure 5 presents a radial cut at r = 4.5 ro for cases 
SOO and SOI. There is a range of angles, ~ 30^ above 
and below the equator, where the MHD wind is slower 
compared to the HD model (SOO). This deceleration is 
a consequence of the magnetic tension that is stronger 
near the equator. At high latitudes, the gradient of both 
thermal and magnetic pressures are responsible for driv- 
ing the wind. Figure 5 also presents the analytical result 
for a non-magnetized wind (dashed line) for comparison 
purposes. The width of the decrease seen in Fig. 5 for 
the SOI curve is dependent on the grid resolution (Opher 
et al. 2004). 

Figure 6 presents the radial cuts for the remaining sim- 
ulations of the first set from ^ = 0^ to 90^. All the curves 
present smaller velocities at the equator {0 = 90^) and 
higher velocities at high latitudes. The difference in ve- 
locities increases as Bq is increased (see Table 2). 



*AWr = Ur{0 = 18°) - Ur{0 = 45°) 



TABLE 3 

Mass-loss rate per unit solid angle of the outflow* 



Name 


mi5(6> = 0°) 


mi5(6/ = 45°) 


S02 


5.1 


5.7 


S03 


6.0 


6.6 


S04 


7.1 


7.5 


S05 


8.2 


8.7 



*mi5(^) = [m(^)/10-i5] Moyr-i 



Another consequence of the latitudinal dependence of 
the Lorentz force, the mass- loss rate per unit solid angle 
of a steady wind 

m{0) =r^p{0)ur{0) (9) 

is also latitude-dependent. According to Eq. (1), rh{0) 
is constant if there is no variation in ^. In our simula- 
tions, this is verified far from the central star. Table 3 
presents the calculated rh{0) for two given co-latitudes of 
the steady-state wind, when it has already reached con- 
stant values. m{0) at the pole is smaller compared to 
other angles, despite the increase of the radial velocity. 
The lower value of m{0) at the pole is due to matter car- 
ried from high latitudes to low latitudes (§3.1.3). As the 
latitudinal increase of density is higher than the corre- 
sponding decrease in the radial velocity, rh{0) is higher 
at low latitudes. 



6 



Vidotto et al. 









t 




7 -3.5 3.5 7 



u| (km/s) 


m 


318 




289 




260 




231 




202 




173 


n 


145 




116 




87 




58 


1 


29 








|u| (km/s) 

326 



65 
33 





f / 


S03] 




1 /-■'■' 
/ /' 








li il /km/c^ 

^ 426 

383 






341 
^ 298 






■ 256 






■ 213 
H 170 




■ 128 

■ 85 












■ '3 

" 









x(ro) 



-3.5 3.5 7 



-3.5 3.5 

x(ro) 







^04^ 






|u| (km/s) 

610 
555 
499 
444 






^ 388 
■ 333 






H 277 

■ 222 

■ 166 












H 55 

™ 


/' / 
^^^^^^^^ 








-3.5 3.5 

x(ro) 

Fig. 4. — Meridional cuts of the steady-state configurations for simulations SOI to SOS. Black streamlines represent the final configuration 
of the magnetic field. Contour plots of the velocity of the fiow are shown in the background. The white line is the Alfven surface. 



220 



200 - 



180 



E 

^ 160 



120 









- - - - Analytical 




soo 




S01 


\ -/* - \ J ✓ - 








1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 





30 60 90 120 150 180 

e(°) 

Fig. 5. — Radial velocity of the wind as function of ^ at r = 4.5 ro, 
for SOO (double-dot-dashed line), SOI (solid line) and analytical 
solution of a non-magnetized wind (dashed line). 



3.1.3. Meridional Flows 




Fig. 6. Same as Fig. 5, but for SOI to SOS, from 6> = 0° to 
= 90°. 



The component of the magnetic force in the direction 
gives rise to meridional flows, bringing matter from both 
hemispheres towards the equator. This causes a density 
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enhancement along the equatorial plane. 

Figure 7 presents the meridional velocity, uq^ for SOI to 
S05 at steady-state, ue increases with Bq, ranging from 
a maximum value of 28 km s~-^ in SOI {Bq = 1 G) to 
339 km s-i in S05 {Bq = 20 G). The last panel of Fig. 7 
presents the flux of matter in the direction, pug^ at 
4.5 tq. It can be seen that as the magnetic field increases, 
there is an increase of the meridional flux of matter. 



TABLE 4 

Location of the neutral point, tn, along the equatorial 

PLANE AS A FUNCTION OF Bq 



Name 


Bo (G) 


rN(ro) 


S02 


5 


3.68 


SOS 


10 


3.90 


S04 


15 


4.01 


S05 


20 


4.12 



3.1.4. The Neutral Point 

The meridional flows are also responsible for compress- 
ing lines with opposite polarities along the equatorial 
plane and due to reconnect ion, a current sheet is formed. 
In our simulations, the reconnection is controlled due to 
numerical resistivity. In order to minimize the effects of 
numerical resistivity, we reflned the grid along the equa- 
torial plane. Ultimately, a more realistic reconnection 
should be included in the global 3D MHD simulations 
such as done by Kuznetsova et al. (2007). The recon- 
nection leads to the emergence of a faster wind in the 
equatorial plane, beyond the reconnection point (Fig. 4). 
Such wind is an artifact of the model and the same be- 
havior was also observed by Roussev et al. (2003). 

The location of the neutral point rjv is given by the 
requirement that the total pressure (magnetic plus ther- 
mal) at the tip of the cusp of closed loops is continuous 



Stt 



(10) 



where p{r^) is the thermal pressure evaluated immedi- 
ately after tn andp(r^) immediately before (PK71). 
Along a magnetic streamline, the magnetic force is null. 
Assuming an isothermal atmosphere, the energy equa- 
tion reduces to 

— + c^ln^ + ^r = ^+^oro, (11) 
2 po 2 

where denotes the isothermal sound speed and the 
use of the subscript "0" indicates that the quantity is 
evaluated at the base of the wind. Equation (10) can be 
written as 



1 



exp ■ 



Cs J 



:,(4)^,(12) 



where we neglected the initial velocity of the flow in com- 
parison with the velocity of the flow at the neutral point, 
7i(rAr), and assumed u{rN) <C Cg. Under these consider- 
ations, Eqs. (10) and (12) lead to 



BivN)' 



:p(r+) 



^ u{rN)^ 



1 



Stt --'-^'^ 2cl "2 
implying that at the neutral point 

u{rN) ^ VAivN) 



p{r^)u{rNf , (13) 



(14) 



where va is the Alfven velocity. I. e., the neutral point 
lies on the Alfven surface. 

Figure 7 presents the Alfven surface (white line) , where 
it can be seen that the above result is confirmed by our 



simulations. In Table 4 the position of the neutral point 
Tat is given for each of the simulations of the first set. The 
cell size at the current sheet is 0.036 ro, so the numerical 
error associated with tat is ±0.018 tq. vn increases as Bq 
is increased, i. e., the "dead-zone" (zone of closed loops) 
is larger for larger ^o- As is increased from 5 G to 
20 G, tn moves from 3.68 ro to 4.12 tq. This is expected, 
since the Alfven surface moves farther from the star as 
Bq is increased. SOI does not present reconnection. 

Inside the closed field line region, particles are trapped. 
This is evidenced in Fig. 8, where we plot the magnetic 
field lines and the vectors of the flow velocity. As can 
be seen, if a given particle emerges inside the closed loop 
zone, it will remain there, because its velocity is not high 
enough to escape the effective potential well created by 
the magnetic tension and gravitational attraction force. 

3.1.5. Energetics 

We would like to estimate the contribution of the dif- 
ferent energies to the total energy content of the system. 

The temperature T of the steady-state coronal wind 
increases with Bq. Figure 9 presents the temperature 
proflles along the polar axis for the cases of the flrst set of 
simulations (in steady-state). For higher Bq^ T increases 
near the base of the wind and it drops slowly for larger 
radii. For S05, the increase in the temperature is max- 
imum among all the cases, reaching c^^ 5.9 MK around 
6 tq. In the solar wind, however, the temperature profile 
is not fiat (Matthaeus et al. 2006), because the bulk of 
the internal energy deposition occurs near the Sun. 

For steady- state, conservation of total energy fiux is 
given by 

V-q = 0, (15) 
where q is the total energy fiux given by 

q = u(^ + ^+p$)+S, (16) 

and $ is the gravitational potential. The terms on the 
right hand side of Eq. (16) are the kinetic, enthalpy, and 
gravitational energy fiuxes, and S is the Poynting fiux 
vector that in ideal MHD is given by 



47r ' 



-u 



47r 



(17) 



Taking a volume defined by a given magnetic fiux tube 
bounded by two cross-sections (see Fig. 10), the net en- 
ergy flux q should be conserved. 

Table 5 presents each parcel of the total energy flux 
that crosses four cross-sectional areas of a magnetic flux 
tube whose central axis is located at ^ ^ 30*^. At 2.2 rg, 
the enthalpy power (i. e., the thermal plus internal en- 
ergy fiuxes integrated over the area of the cross-section of 
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Fig. 8. — Trapping of particles inside the closed loop region for 
S03. Thick lines represent the magnetic field lines and vectors 
represent the fiow velocity, both in the meridional plane. 



the tube at this position) is already the dominant form of 
energy. At 3.0 ro, magnetic power is almost negligible in- 
dicating that the continuous acceleration of the flow from 
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Fig. 10. — Cartoon of the energy flux q along a flux tube. 



TABLE 5 

Percentage amount of energy flux passing 
through four different cross-sections 
along a flux tube for 803 and s05 



Energy Flux 


2.2 ro 


3 ro 


4.2 ro 


6 ro 




S03 


Magnetic 


0.15 


0.07 


0.02 


0.01 


Kinetic 


0.32 


0.59 


1.26 


1.76 


Enthalpy 


97.08 


97.63 


97.50 


97.39 


Grav. 


2.45 


1.72 


1.22 


0.85 


805 


Magnetic 


0.41 


0.15 


0.02 


0.02 


Kinetic 


0.46 


0.85 


1.77 


2.33 


Enthalpy 


98.04 


98.27 


97.72 


97.31 


Grav. 


1.09 


0.72 


0.50 


0.34 



there on is driven by solely the enthalpy. At r = 4.2 ro 
and 6.0 ro, the conversion of enthalpy into kinetic energy 
is responsible for the continuous acceleration of the wind 
until a terminal velocity is reached. 

3.2. The Second Set of Simulations 

The second set of simulations aims to investigate the 
effect of the density on the wind. We vary the initial 
density maintaining the same surface magnetic field ^o 
and surface temperature To adopted in 805. The simu- 
lations we analyze in this section consist of simulations 
S06 to S09 (Table 1). 

We note that as po increases, there is a decrease in the 
wind velocity, as can be seen in Fig. 11, where we plotted 
the radial velocity of the flows along = 45*^. Also, as po 
is increased, the reconnection point moves closer to the 
star. 

3.3. Comparison Between First and Second Sets 




Fig. 
along I 



20 40 60 80 100 

Radial distance along 9=45° 

11. — Radial velocity profiles as a function of the distance 
) = 45° for the second set of simulations. 



Comparing the first and the second set of simula- 
tions, we note that flows with same Po possess the same 
wind velocity profiles (S01/S09, S02/S08, S03/S07, and 
S04/S06). This result can be verified from Eqs. (1) to 
(4). The magnetic field can be written in terms of a 
normalized field B with a normalization constant Bo 

B = B^o 
Doing the the same for the density 

P = PPo- 

In this case, Eqs(l) to (4) are written as 



dt 



puu 



f + V.(pu)=0 



BB 



Po Stt 



Po 



de 
dt 



+V- 




where 



pu^ 



7-1 po 8tt 



(18) 

= PS 
(19) 
(20) 

(21) 
(22) 



Equations (18) and (20) show that the normalized den- 
sity profile {p) and the normalized field profile {B) are 
independent of the values po and Bo adopted at the sur- 
face. Equations (19) and (21) show that the choice of po 
and Bo (or in other words the choice of Po Po/Bq, as 
we assumed the same base temperature) appears explic- 
itly in the momentum and energy equations. According 
to Eqs. (18) to (21), it is expected that flows with same 
Po have the same wind velocity, magnetic field config- 
uration, density profile, etc. Although the normalized 
density and magnetic field profiles are the same if Po is 
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the same in two different simulations, the mass- loss rates 
are different, where the value of the normalization con- 
stant po has to be considered. 

Hence, by analyzing only the velocity of an outflow, 
one cannot equivocally predict its physical characteris- 
tics at its base. Our results indicate that there is a group 
of magnetized flows that would present the same termi- 
nal velocity despite of its thermal and magnetic energy 
densities, as long as /3o is the same. 

This degeneracy can be removed if we consider the 
mass- loss rate. As the velocity profiles for two differ- 
ent magnetized outflows with similar Pq are the same, 
the difference between the mass-loss rates of these winds 
comes from different density structures. To illustrate, 
consider cases SOI and S09 with /3o = 1: both cases 
present the same velocity profile and magnetic field con- 
figuration (although the intensity is not equal); but the 
density, and consequently the mass-loss rate, of S09 is 400 
times larger than SOI in the entire numerical domain. 

3.4. The Third Set of Simulations With 7=1.1 

Throughout the present paper, we have assumed that 
7 = 1.01. However, more realistic magnetized wind mod- 
els compute the heating of the wind caused due to phys- 
ical processes, as is the case of the ID models done by 
Vidotto & Jatenco-Pereira (2006); Falceta-Goncalves et 
al. (2006); Cranmer et al. (2007), among others. How- 
ever, inclusion of such processes, e.g. dissipation of waves 
or turbulence, in a 3D code is very challenging. We thus 
opted to parametrize the energy content of the wind by 
the use of 7. 

We investigate now how the wind structure will change 
if a different 7 is considered. Considering the cases 
S01/S09 and S02/S08, where /3o = 1, 1/25, respectively, 
we performed simulations considering 7 = 1.1. This third 
set of simulations is presented in Table 1 and the merid- 
ional cuts of the steady-state configurations for SO lb and 
S02b are shown in Fig. 12. 

It is worth noting that by adopting a different 7, the ac- 
celeration mechanism of the wind changes, what changes 
the velocity profile of the wind and, consequently, the 
configuration of the magnetic field (compare, for in- 
stance, the steady-state configuration of SOI and SOlb). 
This is expected because, the higher 7 is, the thermal 
acceleration mechanism is less effective (e.g., there is less 
turbulence or dissipation of waves). 

Due to the scaling relations described in §3.3, simula- 
tions SOlb and S09b, with /3o = 1, present the same con- 
figuration of the steady-state wind. The same result is 
achieved for simulations S02b and S08b, with /3o = 1/25. 

4. DISCUSSION AND CONCLUSION 

This is the first study to perform 3D MHD numeri- 
cal simulations analyzing in detail winds of magnetized 
solar-like stars with different /3o. As initial condition, our 
model assumes a bipolar configuration for the magnetic 
field at the stellar surface. This distribution of magnetic 
field, although very simplified, is useful as a first step in 
understanding how the interaction between stellar mag- 
netic field and the wind occur in a 3D self-consistent 
manner. 

In our simulations, we adopt a simpler treatment and 
parametrize the energy content of the wind in terms of 7. 




Fig. 12. — Same as Fig. 4, but for simulations SOlb and S02b, 
where 7 = 1.1. 



This is a first step towards more realistic simulations. It 
should be stressed that if one aims to describe the com- 
plex physics of stellar winds, a more complete description 
of the energy content of the wind has to be taken into 
account. In the solar corona, for instance, it is inferred 
that closer to the Sun, 7 1 (Steinolfson & Hundhausen 
1988), and at 1 AU, 7 1.5 (Totten et al. 1995; Phillips 
et al. 1995). Therefore, 7 is not expected to be constant 
throughout the stellar wind, as we assumed. As the aim 
of our study is to investigate the effects of the magnetic 
fields in the general properties of the stellar wind, the 
inclusion of MHD waves or variable 7 is postponed for 
future work. 

In the first set of simulations where a fixed heating 
parameter 7 = 1.01 was adopted, we showed that for a 
solar-like star, the increase in the magnetic field intensity 
creates faster winds, with general characteristics as the 
PK71 result for the Sun: a creation of a bi-modal stellar 
wind. The final configuration of the magnetic field con- 
sists of a zone of closed magnetic loops at low latitudes. 
On top of the closed loop region lies the neutral point, 
where the reconnection takes place and is the starting 
point of the current sheet that extends for larger dis- 
tances along the equatorial plane. The zone of open field 
lines located at high latitudes of the star fills all the space 
outside the closed loop region and carries along a wind. 
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For a given (3o, the wind is more accelerated in the 
polar regions than in the equatorial region. This differ- 
ence is generated by the latitude- dependent Lorentz force 
which creates a flux of matter directed to the equato- 
rial plane. This generates a density enhancement at low 
latitudes providing a mass-loss rate per unit solid an- 
gle that is also latitude-dependent and increases at low 
latitudes. This meridional flux is intensified for higher 
magnetic fields. The increase in the magnetic field also 
led to hotter winds with higher mass-loss rates. 

We showed that the /^-parameter is a key parameter in 
the structure of the wind in the case 7 is the same for 
all simulations. As /3o was increased back to 1 by chang- 
ing po and maintaining Bq = 20 G in the second set of 
simulations, the wind decelerated and cooled. Compar- 
ing the first and second sets of simulations, we conclude 
that winds with same /3o have i) same velocity profile and 
ii) same magnetic field topology (the neutral point is lo- 
cated at the same position). However, mass- loss rates 
are different as m{0) is dependent on the choice of po. 

By analyzing the effects of a different 7 in the simu- 
lations, we showed in the third set of simulations that 
the heating parameter 7 is, together with the plasma-/^ 
parameter at the coronal base, an important parameter 
in the structure of the wind. However, the steady-state 
configuration of the wind in the third set is different from 
the first and second ones because the thermal force driv- 
ing the wind changes when 7 changes. 

In the simulations we ran, we adopted 7 = 1.01 in the 
first and second sets, and 7 = 1.1 in the third set. Al- 
though 7 is necessary for the thermal acceleration of the 
wind, a wind with a low 7 implies a proportionally high 
Lorentz force, when compared to the thermal force. For 
this reason, the third set of simulations is more asym- 
metrical than the first set (compare SOI and SOlb, S02 
and S02b), since the Lorentz force is ^-dependent. An- 
other aspect is that in a purely HD coronal wind, the 
thermal force is the main responsible for the accelera- 
tion of the wind. This implies that by adopting a low 
7, the wind has less internal energy available, being less 
thermally accelerated. In the cases SOI (7 = 1.01) and 
SOlb (7 = 1.1), the same situation holds: case SOlb 
is less accelerated than case SOI. However, cases S02 
(7 = 1.01) and S02b (7 = 1.1) show a different behavior. 
The present simulations indicate that a change in gamma 
from 1.01 to 1.1 decreases the terminal speed by 43% in 
cases SOI /SOlb and increases 8% in cases S02/S02b. This 
different behavior is because in case S02b, the Lorentz 
force provides a significant acceleration mechanism for 
the wind, i.e., the Lorentz force is proportionally more 
important than the thermal force. This illustrates the 
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need for perform full 3D MHD simulations in order to 
assess the importance of (3 and 7. 

The effective temperature of solar-like stars is orders of 
magnitude lower than the coronal temperature adopted 
in our models. In the present work, we did not take into 
account what is heating the wind from photospheric tem- 
peratures up to coronal temperatures. Instead, we placed 
the inner boundary of our model at the coronal base. 
For a more realistic study, however, the energy equation 
should be solved starting from lower layers. Hence, the 
plasma- /5 parameter at the coronal base would then be 
incorporated as resultant from the energy interactions 
happening in these lower layers rather than being a free 
parameter at the coronal base, as in our models. Mod- 
els and ID simulations that resolves the energy balance 
in the chromosphere-corona transition region have been 
calculated recently for the solar wind (Suzuki & Inut- 
suka 2005; Cranmer et al. 2007) and stellar winds (Suzuki 
2007; Cranmer 2008). If incorporated in 3D global sim- 
ulations, these models can provide physical insights on 
the processes that accelerate and heat stellar winds. 

The present work considered moderate f3o. In the pres- 
ence of a weak magnetic field (high /^o), the wind is en- 
ergetic enough to drag the field lines with it, leading 
to a configuration of radial magnetic field lines. Conse- 
quently, a spherical expansion of the corona is expected, 
similar to a purely HD case. 

Here, we deal with a class of non-rotating stars. The 
validity of our results should be tested also in rotating 
objects, since a more generalized description of magne- 
tized outflows from solar-like stars requires the inclusion 
of rotation. The detailed interaction including rotation 
requires detailed 3D MHD simulations. A future work 
will extend the present work to this direction. 
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